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1 Introduction 


The theory of differential invariants nowadays undergoes active development and is widely 
used for integration in quadratures and for order lowering of ordinary differential equations, 
and also for description of classes of invariant differential equations |^, . In the theory of 

differential invariants a major role is played by various versions of the conjecture on the hnite 
basis of differential invariants that could be non-rigorously formulated in the following way: 
for an arbitrary group G of local transformations there exists such finite set of differential 
invariants that every differential invariant of the group G can be represented as a function of 
these invariants and their derivatives. A statement of such type (for one-parameter group of 
local transformations in the space of two variables) was proved by S. Lie himself (see also 
Bi) and soon afterwards it was essentially generalized by A. Tresse . The recent progress 
in this direction is due to the works by L.V. Ovsyannikov 0 and P. Olver B, I, I, [Tod , where 


the notions of the operator of invariant differentiation, differential invariant coframe etc. are 
introduced, and results on rank stabilization of the prolonged group action and on the 
estimates for the number of differential invariants are obtained. There is also considerable 
number of papers devoted to the search for the differential invariants of specihc groups (see 

e.g. 0, O, H, [ 


li)- 


In the present paper we study the differential invariants for one-parameter group of local 
transformations in the space of n independent and m dependent variables (m, n G N). (The 
case n = m = 1 was considered in |I6| , 0] and in our work |0]. The results for the case 
n = 1 with arbitrary m G N were published in |0].) The importance of this problem 


stems from its being a part of the problem of search for differential invariants of a group of 
arbitrary dimension i§- Our generalization of the Lie theorem on differential invariants 
of a one-parameter group of local transformations was done not only with respect to the 
number of independent and dependent variables, but also with respect to strengthening of 
its statements. It was proved that if a differential invariant is known, then a complete set of 
functionally independent differential invariants can be constructed through one quadrature 
and differentiations. As a side product, we also present some results on the existence of 
basis of differential invariants rational in higher jet coordinates. A link between differential 
invariants and integration of systems of Riccati-type equations was analyzed. Let us note 
that results of this paper can be generalized for some classes of multiparameter groups of 
local transformations (or Lie algebras of differential operators). 


2 Generalization of Lie theorem 
on differential invariants 

Let Q = f°‘{x,u)dx^ + r]^{x,u)dui be an inhnitesimal operator of a one-parameter group G 
of local transformations which act on the set M C J(o) = X x U, where X ~ M"’ is the 
space of independent variables x = {xi,X 2 , ■ ■ ■ ,Xn) and U ~ is the space of dependent 
variables u = ... ,m"^), is a prolongation of the action of the group G for the 

subset M{r) = M X x x ■ ■ ■ x of the jet space J{r) = X x U(r) of r-th order 
jets over the space X x U (here U(^r) = U x x x • • • x r ^ 1, is the r-th 
prolongation of Q B; ^). A function J: M^^.) —> M is called a differential invariant of the 
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order r for the group G (or for the operator Q) if it is an invariant of the prolonged action 
of (of A necessary and sufficient condition for the function I to be an r-th order 

differential invariant of the group G is the equality = 0. 

Here and below, if not otherwise stated, the indices a, b, c, d run from 1 to n, indices i, 
j, k, I run from 1 to m. The summation over the repeated indices is understood. 

Let / = I{x,u) = P{x,u),..., be a complete set of function¬ 

ally independent invariants (or a universal invariant 0) for an operator Q, and J{x,u) 
is a particular solution of the equation QJ = 1. Then the functions I^{x,u), P{x,u), 
..., u) and J{x,u) are functionally independent. Let us make a local change 

of variables: Hc = P{x,u), c = l,n — 1, yn = J{x,u) are new independent variables, and 
yi _ p+^-^(^x, u) are new dependent variables. In terms of variables y = {yi,y 2 ,, yn) and 
V = ... ,n™') the operator Q has the form dy^. Thus for any r ^ 1 the form of the 

prolonged operator coincides with Q = dy^, and therefore 


y = (2/1,2/2, • • -,2/^-1), 

Q\a-\p 


Up) — 


dyTdyT • • • ^ 2 /^ 


aa£hlU{0}, |a| =y~~^ Oia^r 

a=l 


(here = n*, if |q:| = 0 ) form a complete set of its functionally independent invariants, 
and {y,V(^r)) is a universal invariant of the group G^^\ (Functional independence of the 
components y and V(^r) is obvious, as {y,V(^r)) is a set of variables in the space J(r)-) This 
means that {y,v) is a fundamental set of differential invariants for the operator Q, i.e. any 
differential invariant of the operator Q can be represented as a function of y and v and of 
the derivatives of v with respect to operators of G-invariant differentiation. These operator 
coincide here with the operators Dy^ = dy^ -|- Vy^d^^i + Vy^y^dviy^ H— ■ of total derivatives with 
respect to the variables 2 /a- 

Let US go back to the variables x, u. In terms of these variables 


D 


Vc 


D 


yn 


{-IY+- D{I^,d=T^i,d^c, J) ^ 

A Dl^Xh, 6=1,71, b^a) 

(_l)n+a ^ 

A D(^Xf), 6=1,71, b^a) 


c = 1 , n — 1 , 


( 1 ) 


where = dx^ + u^dni + ■ is the operator of total derivative with respect to 

the variable Xa, and 


D{I'^,d=l,n-l,d^c, J) D{I‘^,d=l,n-l) ^ d=l,n-l, J) 

D(^Xh, 6=1,71, b^a) D(^Xfy, 6=1,71, b^a) Dl^Xb, 6=1, 71 ) 

denote Jacobians (of total derivatives) 

of the functions d = 1 , n — 1 , d 7 ^ c, J with respect to the variables 6 = 1 , n, 6 7 ^ a, 
of the functions d = l,n — 1 with respect to the variables Xb, 6 = 1 ,n, 6 7 ^ a, 
of the functions /‘^,d = l,n — 1, J with respect to the variables Xb, b = 1, n, 
respectively. 

As a result, we arrive at the following theorem. 
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Theorem 1. Let I{x,u) = P{x,u),..., ^{x,u)) be a universal invariant of 

an operator Q and J{x,u) be a particular solution of the equation QJ = 1. Then functions 




r)oii n«2 

yi y2 ■ ■ ■ yn ^ 


- 1 / 


where c = — G M U {0}, X]a=i ^ operators Dy^ are determined by 

the formulae w, form a complete set of functionally independent r-th order differential 
invariants (or a universal differential invariant) for the operator Q. 

Corollary 1. For any operator Q there exists a complete set of functionally independent r-th 
order differential invariants, where every invariant is a rational function of the variables 
(a = (oi, 02, • • • , C(n), Oa G N U {0}, 0 < X]a=i ^ of the jet space J(^r) with coefficients 
depending on Xa and . 

Corollary 2. If I = {I^{x,u), P{x,u),..., is a universal invariant for the 

operator Q and J = J{x,u) is a particular solution for the equation QJ = 1, then the 
functions 


ri+n—l 

^y-n^ 


d=l,n-l, d^c, J) ^ ji+n-l 

A b=l,n, 

A D{xi)^ 6=1,n, 


( 2 ) 


(c = l,n — 1) form a complete set of functionally independent differential invariants having 
exactly order one for the operator Q. 

Note that if a universal invariant I for the operator Q is known, then a particular solution 
for the equation QJ = 1 may be easily found via one quadrature. E.g. if for some hxed a 
ff 7 ^ 0 , we have a particular solution 


J{X, U) = J dXa/^'^iX^ . . . , X--\Xa, , X^, U\ . . . , U^), 

where Xb = X^{xa,C), b ^ a, = U^{x,C) is the solution for the system of algebraic 
equations I{x,u) = C ■.= {Ci,C 2 , ■ ■ ■ ,Cm+n-i) with respect to the variables xt, b ^ a, uf 
and after the integration it is necessary to perform an inverse substitution C = I{x, u) (there 
is no summation with respect to a in this case). Likewise, when 77 * 7 ^ 0 for some fixed i, then 
we can assume 


J{x, u) = f duQp^Xf ...,X^,U\..., U^-\ u\ f/*+i, 

(there is no summation with respect to i), where Xb = X^(m*,C'), = U^{u\C), j 7 ^ i, is 

the solution of the system of algebraic equations / {x, u) = C with respect to the variables 
Xb, ufj^ i. 

Thus, the following theorem holds. 

Theorem 2. If a universal invariant is found for the operator Q, then a complete set of 
functionally independent differential invariants of arbitrary order may be constructed through 
one quadrature and differentiations. 
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3 Invariant Differentials 


Let us introduce a notion of an invariant differential that is a particular case of a more 
general notion of a first-order contact invariant differential form in the jet space [Q. 

Definition. A differential dW{x,u) will be called invariant with respect to a group G (with 
operator Q), if it does not change under action of transformations from the group G. 

A criterion for invariance of a differential is an equality dQW{x,u) = 0. Two essentially 
different cases are possible: 

1) the function W{x,u) is an invariant of the operator Q, i.e. QW{x,u) = 0; than the 
differential dW{x,u) is automatically invariant with respect to the operator Q {invariant 
differential of the first type)] 

2) the function W{x, u) is not invariant under Q, while the differential dW (x, u) is {invariant 
differential of the second type)] then QW{x,u) is a non-zero constant. 

If a set of functions I{x,u) = {I'^{x, u))g^ i jn+n-i ^-nd J{x,u), that determine a universal 
invariant of the operator Q and the invariant differential of the second type is type, then all 
such sets may be found according to the formulae 

I{x,u) = F{I{x,u)), J{x,u) = J{x,u) + H{I{x,u)), (3) 

where F = {F^, F"^,..., F"^~^^~^) and H are differentiable functions of their arguments, 
\dF/dI\ 7 ^ 0. The formulae @ determine the equivalence relation on the set Ai of collec¬ 
tions of m -|- n smooth functions of m -|- n variables with a non-zero Jacobian. We denote 
the corresponding set of equivalence classes as A4/Q. 

Proposition. There is a one-to-one correspondence between Ai/Q and the set of non-zero 
operators {Q} in the space of variables {x,u): the set {{I{x,u)] J{x,u))} of solutions of 
the system QF = 0, q = l^m + n — 1, QJ = 1, where F are functionally independent, 
is an element of the set Ai/Vt, and vice versa, if {I{x,u)] J{x,u)) is a representative of 
an equivalence class from Ai/fl, then the system QF = 0, q = l,m F n — 1, QJ = 1 
is a determined system of linear algebraic equations with respect to the coefficients of the 
corresponding operator Q. 


4 The case n = 1 

Let us consider in more detail the case of one independent variable x {n = 1), for which it is 
possible to obtain a more compact formulation of Theorem 1 and of its corollaries, and also 
to obtain some additional results. 

Theorem 1'. Let I = I{x,u) = {F{x,u), F{x,u),..., I"^{x,u)) be a universal invariant 
of the operator Q and J{x,u) be a particular solution of the equation QJ = 1. Then the 
function 

F{x,u), s = T~f, 
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where = dx + u^d^i + u^xx^ui. + • ■ ■ is the operator of total derivative with respect to the 
variable x, form a complete set of r-th order functionally independent differential invariants 
(or a universal differential invariant) of the operator Q. 

Corollary 1'. For any operator Q there exists a complete set of n-th order functionally in¬ 
dependent differential invariants, where every invariant is a rational function of the variables 
ul, . .., of the jet space with coefficients depending on x and u\ 

Corollary 2'. If I = (/^(x, u), P{x, u),..., /™(x, u)) is a universal invariant of the operator 
Q and J = J{x,u) is a particular solution of the equation QJ = 1, then the functions 


^(1) = ^(i)(^>“(i)) 


dP 

dT 


DxP 

DxJ 


Px + ijiUj 
Jx + Jui''^x 


( 4 ) 


form a complete set of functionally independent differential invariants of exactly first order 
for the operator Q. 

Corollary 3. The components of universal differential invariants having exactly order one 
of the operator Q may be sought for in the form of fractional-linear functions of the variables 
ul, of the jet space with coefficients depending on x and u\ 

Corollary 2' may be restated using the notion of the invariant differential. 

Corollary 4. The ratio of invariant differentials of the operator Q of first and second type 
is its differential invariant of exactly first order. If dP, dP, ..., dl"^ form a complete 
set of independent invariant differentials of the first type for the operator Q, then its ratio 
with its invariant differential of the second type exhaust functionally independent differential 
invariants of exactly first order for the operator Q. 

Corollary 5 (Lie Theorem) [^, |^, Let n = m = 1, I{x,u) and fiifix,u,Ux) are 
differential invariants of zero and of exactly first order for the operator Q. Then the functions 


P 


hi)) 


dP 



s = 1, r — 1, 


form a complete set of n-th order functionally independent differential invariants for the 
operator Q. 

The operators of G-invariant differentiation for the case of one independent variable are 
traditionally sought for in the form 


V 


1 


Dx, 


where P is a differential invariant for the group G (see e.g. Corollary 5). Therefore, for the 
construction of an arbitrary-order universal differential invariant for a one-parameter group 
of local transformations by means of the operator of G-invariant differentiation of this form 
it is necessary to know m -|-1 functionally independent differential invariants for the group G 
of the possibly minimal, or m functionally independent zero-order differential invariants (or 
simply invariants), and one exactly hrst-order differential invariant. The algorithm suggested 
in the Theorem 1' allows to avoid direct construction of differential invariants. 
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Example 1. (Cf. ^ .) Let n = m = 1 and G = S0(2) be a group of rotations acting on 
X X U ~ with an infinitesimal operator Q = udx — xdu- I = \/+ v? is an invariant 
of the group G (of the operator Q), whence (in notation from the proof of Theorem 2) 
U{x, G) = ±VC'2 -xT Then 


J = ± 


dx 


VC2 - 


, . a; , . X 

= ± arcsm — = ± arcsm 


G 


yj x"^ + 




(here we put the integration constant to be zero), whence 


Ad 


d-X T lu^X 

Ae T Ju^X 


X + UUx 
—U + XUx 


\/a;^ + 


or J(i) 


X + UUx 
—u + XUx 


is a first-order differential invariant for the operator Q. 


5 The Standard Approach 

and Integration of Riccati-Type Systems 


Within the framework of direct method the differential invariants having exactly first order 
are found as invariants of the first prolongation 


+ (he + vtx - C 


Uu 




)du 


of the operator Q, or as the first integrals of the corresponding characteristic system of 
ordinary differential equations 

dxa du^ du^ 

vc+vtx - ex A - ’ 

that depend not only on x and u, but also on the other the variables of the space J(i). (Here 
is a variable of the jet space J(i), which corresponds to the derivative du'^/dxa] lower 
indices of functions stand for the derivatives with respect to the corresponding variables; 
there is no summation over a, c, i and k in the latter equation). Integration of the system 
(^) is, as a rule, a highly cumbersome task. If a universal invariant I{x,u) for the operator 
Q is known, then it amounts to integration of Riccati-type systems of the form 


du’^ 


dx 


a 


Gb r^k 

-^uP,ul + YuP, 

Ca ^ ^ ta ^ 


_L Aa:c 


U=U(Xa,C) 

xa=X‘^{xa,C), 


if 7 ^ 0 for some fixed a, or 


( 6 ) 


du^ 

du^ 


r^i 


uiul + ^ul 


e 

V 




x=X{u\C), , 

u’-=U\u\C), 


( 7 ) 


if 77 * 7 ^ 0 for some fixed i. Here Xd = X‘^{xa,G), d ^ a, u = U{x,G) and x = X{u\G), 
u^ = U\u\G), I 7 ^ i, are solutions of the system of algebraic equations I{x,u) = G with 
respect to the variables Xd, d ^ a, u and x, u\ I 7 ^ z, respectively. The constants G = 


7 


















(Cl, ( 72 ,..., Cm+n-i) in the systems and (]^) are considered as parameters. The case 
p* 7 ^ 0 conld be rednced to the case 7 ^ 0 by means of the locns transformation: 


Xa = u\ 


Ua = —, 

ui 


Xd = Xd, n* = Xa, =u\ a, i, 


< 

~i < 

~l 1 '^d 1 

i ’ 

«a = —, 

Ud = Ud- —Ua 

K 

K 

K 


For this reason we will consider in detail only the case 7 ^ 0. 

The method we snggested in Corollary 2 for hnding differential invariants having exactly 
hrst order, nnlike the standard method, allows to avoid direct integration of systems of 
Riccati eqnations or ([^ and to hnd a solntion throngh one qnadratnre and differentiation. 
This resnlt means that in the case of known universal invariant I{x, u) for the operator Q the 
systems and are always integrable by means of one quadrature. Really, the general 

solntion for the system (j^) conld be given explicitly by m non-linked systems of linear 
algebraic eqnations 




U=U{Xa,C) 
xa=X<^{xa,C), d^a 


0 , 


where P = P~^'^ ^ Cjd/'^ + CjnJ, Cn, are arbitrary constants. To write the solntion 

in the explicit form, we introdnce some additional notation: 


X = {Xd)d=l,d^a^ ^ ^ = Xa, 

jx ^ ju ^ 

/^x _ ! \n—1 /~iu _ !\m 

O' — y^d)d=ly ^ — [^j+n-l)j=l, 

a = (c,,)™, , c" = {c,a)fLi , i = p + c'p + c"j. 

Then the general solntion for the system is given by the formnlae 


^l)T=i h=i 


n 



U=U{Xa,C) , 
X=X{Xa,C) 


or 


«)-1 = U,- Uc-^X-lX, + H{{C' + C"Jc-^)X-lX, - C"Jc^), 
iLi,,^, = Uc-^X-l - H{C' + C"Jc^)X-l, 


where H = {Uc^ - Uc^X-^Xcu){E + C"Jc- - (C' + C"Jc^)X-^Xcu)-\ E is the m x m 
nnit matrix; the signs of vector-fnnctions with lower indices of the sets of variables designate 
the corresponding Jacobi matrices. To ensnre the existence in some neighborhood of a hxed 
point u^) of all inverse matrices, that are mentioned above, it is snfficient to consider 
constants Cn, to be small and to perform a (preliminarily determined as non-degenerate) 
linear change in the set of invariants for the matrix I(^x,u){x^,u^) to be the nnit matrix. 

If we pnt Cib = 0, then we obtain a particnlar solntion 


= C, - Uc-^X-lX,, 


K)T=lUb^a = Uc^X-l 




The solution of the system in explicit form should be written separately for the case n = 
1. As in this case u = U{x, C) is the general solution of the system du^/dx = rj^{x, u)/^{x, u), 
then it is easy to verify that Ux = Ux{x, C) is a particular solution for the system (||) (here, 
as in (^, (7 is a set of parameters). The general solution for the system (H) has the form 


'll 2 


= -{Iu-C®J^)-\lx + CJx) 


u=U{x,C) 


= U,- UciE + C® JcY'CJ, 


( 8 ) 


(in the latter equality the substitution x = z, u = U{z,C) was performed), where E is 
the unit matrix of the dimensions m x m, C = (Ci, C 2 , ■ ■ ■, CmY is a column of arbitrary 
constants, A = (/*,), 4 = (/*), Y = Uc = C®Ju = C®Jc = i&Jci). 

The inverse matrices in (^ always exist for sufficiently small Ci. 

Example 2. Let n = m = 1, Q = exp(— a; — u){dx + udu)- I{x,u) = uexp{—x) is an 
invariant for the operator Q, whence U{x,C) = C'exp(a;). Then 


J 


dx 

exp(—X — C exp(x)) 


— exp (C exp (x)) 

o 


exp(x + u) 
u 


(here we set the integration constant to be zero). Therefore, 


Y) = 


Ix + luUx _ exp(-2x - u)u^{ux - u) 

Jx ~\~ Ju^X ^ tltlx — 


or J(i) = 


Ux — u 


U + UUx — u. 


■ exp(—u) 


is a first-order differential invariant for the operator Q. System (^) for the operator Q 
consists of one equation which has the form 


—= Ux + {2 — C exp(x))-Ua; — C exp(x). 
dx 

The function v = C exp(x) is its particular solution. The general solution of this Riccati 
equation is given by the formula 


Ux = C exp(x) 


C2exp(2x) 

C exp(x) — 1 -|- Cexp(—(7exp(x)) 


where C is an arbitrary constant. 

Example 3. Let n = m = 1, Q = xu{xdx + kudu), G M. I{x,u) = ux~^ is an invariant 

^ ' I ’It 


if k = —1, 
-, if k ^ —1 


for the operator Q, whence f/(x, C) 


" In X In X 

f dx 

C XU 


(/c+l) 


i {k + l)C 


(here we set the integration constant to be zero). The corresponding Riccati equation has 
the form 


dUx 

dx 


J_ 2{k-l) 

Cxt ” X 


Ux + kCx 


k-2 
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The function Ux = kCx^ ^ is its particular solution. The general solution of this Riccati 
equation is given by the formulae 


Ux = -^ I 1 + 

X‘ 


Ux = Cx 


k-l 


C — \nx 
k + 1 


1 + Cx^^^ 


lik = — 1 , 

, if fc 7^-1, 


or 


where C is an arbitrary constant. 


Remark. For well-known transformation groups on the plane (i.e. n = m = 1) integrability 
in quadratures of equations (|^) and (|^) as a rule obviously follows from the form of these 
equations. For instance, when = 0 or = 0, they are a linear equation or a Bernoulli 
equation respectively. If G is a one-parameter group of conformal transformations, then 
ix = Vu and = —Vx, and therefore in equations (|^) and (|^ variables are separated: 


dv 
-I- 1 



u=U{x,C) 


and 


dv 

v‘^ + 1 


'IX , 

—du 
V 


x=X(u,C) 


Example 4. Let n = 1, m = 2, Q = exp(—a; — u^— u^){dx -l- u^du^ + u^du'^). /^(x, u^) = 

M^exp(—x) and P{x,u^,u'^) = ■u^exp(—x) are invariants for the operator Q, whence 
U^{x,C^,C‘^) = G^exp(x) and f/^(x,G^,G^) =G^exp(x). Then 

exp((G^ + G2)exp(x)) 


j{x,C\C^) = 


dx 




exp (—X — {C^ + G^) exp(x)) 

(here we set the integration constant to be zero). The corresponding Riccati-type system 
has the form 
du^ 

-^ = {ul + ul) ul+{2- exp(x)) ul - exp(x)n^ - G^ exp(x), 

J -.2 

= («x + «x) exp(x)M^ -t- (2 - G^ exp(x)) ul - exp(x). 

, the general solution of this system is given by the formula 



= exp(x) 


G^ 

G^ 


+ 


(Gi + G2)exp(2x)J(x,Gi,G2) 


1 - (Gi + G2) (exp(x) - (Gi + G2)"^) J (x, C\ C^) 



where G^, G^ are arbitrary constants. 

Example 5. Let n = 1, m = 2, Q = exp(M^ -|- u‘^){dx + u^du^ — u^du-^). Then, 
(x,G^,G^) = G^ cosX -|- G^ sinX, 

(x, G^, G^) = —G^ sin X -|- G^ cos x, 

j(x, G^,G^)= f exp (—(G^-I-G^) cosx — (G^ — G^) sinx) dx. 
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The corresponding Riccati-type system has the form 


dul 

dx 

dul 

dx 


+ ul) ul + (—sin X + C'^ cos x) ul + (—sin x + cos x + l) ul, 

{ul + ul) ul — {C^ cos X + sin X + l) ul — (C^ cosx + sinx) ul- 


It follows from (j^) that the general solution of this system is given by the formula 

/ “x A _ / —C^ sin X + cos \ , 

I ^2 ) cosx — sinx j 

exp (-(C^ + (T^) cosx - ((7^ - C^) sinx) / cos x + sin x 

1 - C^Jci - C^Jc2 I sin X + cos X 


where C^, are an arbitrary constant. 
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